Micro canonical entropy for small magnetisations 



Hans Behringer 

Institut fiir Thcoretische Physik I, Univcrsitat Erlangcn-Niirnbcrg, D - 91058 
Erlangcn. Germany 

E-mail: Hans . BehringerOphysik . uni-erlangen . de 

Abstract. Physical quantities obtained from the microcanonical entropy surfaces 
of classical spin systems show typical features of phase transitions already in finite 
systems. It is demonstrated that the singular behaviour of the microcanonically 
defined order parameter and susceptibility can be understood from a Taylor expansion 
of the entropy surface. The general form of the expansion is determined from the 
symmetry properties of the microcanonical entropy function with respect to the order 
parameter. The general findings are investigated for the four-state vector Potts model 
as an example of a classical spin system. 
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1. Introduction 

The macroscopic behaviour of a physical system in thermodynamic equihbrium is related 
to the microscopic properties by statistical mechanics. The basic quantity in this 
connection is the density of states depending on the macroscopic quantities of interest. 
For a classical spin system as a model system to describe magnetic properties these 
macroscopic variables are the energy and the magnetisation. In the traditional approach 
to the statistical description of phase transitions the density of states is Laplace- 
transformed to the partition function which is a concept of the canonical ensemble 
and the physics is deduced from the corresponding potential, the Gibbs free energy. In 
recent years an alternative approach to analyse phase transitions within a statistical 
framework has been developed in the microcanonical ensemble (Gross, 1986a; 1986b; 
2001; Hiiller, 1994; Promberger/Hiiller, 1995; Gross et al., 1996). Apart from works 
about general questions of the equivalence and inequivalence of the various ensembles 
(Lewis et al., 1994; Dauxois et al, 2000; Barre et al., 2001; Ispolatov/Cohen, 2001) 
and of the thermodynamic limit in the microcanonical ensemble ((Kastner 2002) to 
cite a recent work) also second order phase transitions have been studied recently 
(Kastner et al., 2000; Hiiller/Pleimling, 2002). Ways to extract critical exponents from 
microcanonical quantities have been suggested (Kastner et al., 2000; Hiiller/Pleimling, 
2002). In these works is was numerically demonstrated that the microcanonically defined 
physical quantities like the microcanonical order parameter or the susceptibility exhibit 
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singularities already in finite systems. However, these numerical studies suggest that the 
microcanonical singularities are characterised by classical critical exponents in contrast 
to the non-trivial exponents showing up in the thermodynamic limit. 

The aim of the present paper is to provide an analytic understanding of the 
characterisation of the singularities in physical quantities of finite microcanonical 
systems by classical exponents. To this end the entropy surface is expanded into a Taylor 
series in terms of its natural variables energy and magnetisation and the asymptotic 
behaviour of the macroscopic quantities is investigated. This procedure is similar to the 
Landau expansion of the free energy of the infinite system although a striking difference 
has to be stressed. The Landau approach is an approximation of the free energy of the 
infinite canonical system whereas the treatment presented below is exact for the finite 
microcanonical system. 

The work is organised in the following way. The second section contains a short 
survey of the microcanonical approach to statistical mechanics. In particular the 
basic concepts and findings of the investigation of second order phase transitions are 
summarised. In the third section the appearance of the classical critical exponents is 
related to the analyticity of the thermodynamic potential of finite systems. The fourth 
section contains general considerations about the Taylor expansion of entropy surfaces 
with a Cnv order parameter symmetry. In addition an investigation of the entropy 
surface of finite vector Potts models with four states is presented. The analysed data 
are obtained from Monte Carlo simulations. 

2. Thermostatics in microcanonical systems 

The analysis of statistical properties of classical spin systems within the microcanonical 
ensemble starts from the density of states 



Here it is assumed that the system is defined on a d-dimensional hypercube of extension 
L. The Hamiltonian Ti gives the total energy of the microstate a and the magnetisation 
operator A4 its total magnetisation. In general the magnetisation is a multicomponent 
object M = (Ml, . . . , M„). The microcanonical entropy is obtained by taking the 
logarithm of the density of states: 



To compare different system sizes the extensive factor L'^ is divided out to give the 
specific energy e = E/L'^ and the specific magnetisation m = M/L'^. The dependence 
of the physical quantities on the system size L is suppressed in the following. The 
statistical properties of an isolated physical system are deduced from the entropy surface. 
The microcanonical spontaneous magnetisation of a finite system for a fixed energy is 
defined by 
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hence msp(e) corresponds to the magnetisation with the maximum entropy for a given 
energy. 

The non-vanishing multicomponent spontaneous magnetisation vector defines a 
direction in the order parameter space: 

msp(e) = (m(p)(e), . . . ,m(p)(e)) = |msp(e)|mo (4) 

where tuq is a fixed unit vector with several possible orientations below Cc- Note however 
that it might be possible that mo and consequently the associated orientations are 
themselves functions of the energy below Cc- Let Go be the symmetry group of the 
microcanonical entropy with respect to the magnetisation components, i.e. 

s{e,g{m)) = s{e,m) (5) 

for all transformations g in Go- If the spontaneous magnetisation rusp is finite the 
symmetry group Gq of the entropy is broken down to the subgroup G < Gq that leaves 
the direction mo invariant (Behringer, 2003). 

In the microcanonical ensemble of isolated systems no external magnetic field 
appears. However, magnetic fields can be defined as the conjugated variable of the 
magnetisation components giving rise to the relation^ 

d d 
/ii(e,msp(e))— s(e,msp(e)) = -^^s(e, msp(e)) (6) 

for the zth component in equilibrium. 

The susceptibility of the system is related to the curvature of the entropy along 
the magnetisation direction. For a multicomponent magnetisation the curvature at the 
spontaneous magnetisation is related to the Hessian 

Hki{e,m^p{e)) = — s(e,msp(e)) (7) 

omkOmi 

of the entropy and the susceptibility is consequently a tensor which is given by 

x!?(e) = -^.(e,m.p(e)) (//"^(e, m.p(e)))^^. . (8) 

From a geometrical point of view the specific heat of the system is connected to 
the curvature along the energy direction. At the spontaneous magnetisation the specific 
heat is given by 



,io)^e) = -^j0^^^^^^^. (9) 
^s(e,msp(e)) 
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Alternatively one can consider a specific heat that is more directly related to the 
canonical viewpoint. There the specific heat is calculated from the canonical entropy 

5can = -^^?(T,/l) (10) 

X The sloppy but more convenient notation ■§^f{xo) for -§^f{x)\x=XQ is used in the following, upper 
indices are not discriminated from lower ones. 
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obtained from the Gibbs free energy g{T, h) of the system in equihbrium as 

d 

Defining the equihbrium entropy by 

5(e) = s(e,msp(e)) (12) 
one obtains the alternative specific heat 

c(e) = . (13) 

The specific heat c is different from c^^^ as nigp is first plugged into s and the derivative 
with respect to the energy is worked out afterwards. In contrast, c^^^ is evaluated 
according to footnote \. 

From a statistical point of view a phase transition from a disordered high-symmetric 
macrostate to an ordered low-symmetric macrostate is defined to take place at a non- 
analytic point of the corresponding thermostatic potential describing the properties 
of the physical system. Yang and Lee showed that the grand-canonical thermostatic 
potential, i.e. the Gibbs free energy, is an analytic function for all finite system sizes 
(Yang/Lee, 1952; Lee/Yang, 1952). This means that a phase transition can only occur 
in the thermodynamic limit of an infinite system, according to the definition above phase 
transitions are not possible in finite systems. 

However the microcanonical equilibrium quantities exhibit typical features of 
phase transitions (Kastner et al., 2000). For instance the microcanonical spontaneous 
magnetisation © displays the characteristics of spontaneous symmetry breaking. This 
spontaneous magnetisation defines the order parameter. For the high energy phase the 
microcanonical order parameter is zero reflecting a phase with high symmetry. The 
abrupt emerging of a finite order parameter for low energies indicates the transition 
to an ordered phase with lower symmetry. The energy at which this onset occurs 
defines unambiguously the critical energy Cc of the finite microcanonical system. The 
susceptibility (jH)) of the system diverges at this energy. 

The use of the expressions critical and phase transition to describe the behaviour of 
finite microcanonical systems is somewhat problematic as they are commonly used for 
certain properties of the infinite system. Nevertheless the appearance of features of the 
microcanonical quantities which are also found in infinite systems at the transition point 
suggests its usage. Another delicate point is the discreteness of the physical quantities 
of discrete spin systems. There the used language refers to continuous functions that 
describe the discrete data most suitably. In the microcanonical analysis of continuous 
spin systems as e.g. the xy (Richter et al., 2004) or the Heisenberg model this concern 
does not exist. 

The existence of a nontrivial spontaneous magnetisation is in general related to 
the appearance of a convex dip in the microcanonical entropy. This is different in the 
canonical ensemble where the curvature along the natural variables is connected to the 
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mean square deviation of a physical quantity. The curvature has therefore a well-defined 
sign. Although the physical quantities show singularities at the transition point that are 
typical of phase transitions there is no phase transition as long as the microcanonical 
system is finite. The microcanonical entropy of finite systems is an analytic function 
of its natural variables and hence a phase transition in the above defined mathematical 
sense does not take place. 

The singular behaviour of the order parameter with respect to the energy can, 
however, be described by means of a critical exponent (Kastner et al., 2000). For a 
physical quantity Q{e) of the finite system this critical exponent, say k, is defined by 



The microcanonical critical exponent (3 of the spontaneous magnetisation in finite 
systems turns out to have the classical value 1/2, the corresponding exponent 7 of 
the susceptibility is 1 for all system sizes. Both definitions and (fT!^ of the specific 
heat are characterised by the exponent a = 0. Whereas the specific heat c^^^ has a 
kink at the critical energy, c exhibits a jump at Cc with different values c(e^) 7^ c(e~). 
The kink of c*-°^ at the critical energy with c*^°''(e^) = c^^\e~), but different left-handed 
and right-handed derivatives at Cc, however, is not a cusp-singularity according to the 
definition of (Stanley, 1972) as the higher order derivatives of c^^^ do not diverge at Cc- 
Note also that both quantities tend to the same function in the thermodynamic limit. 

3. Taylor expansion of the microcanonical entropy 

For a finite system with the typical characteristics of a continuous phase transition 
the microcanonically defined order parameter near the critical point vanishes like a 
square root. This behaviour can be understood in a Landau type of approach to the 
description of phase transitions in the microcanonical ensemble. To this end the entropy 
is expanded into a Taylor series as a function of the order parameter components. A 
proof of the analyticity of the microcanonical entropy is still lacking. Nevertheless in 
view of the analyticity of the canonical thermostatic potential it seems obvious to use this 
assumption. Note that this approach is exact provided the entropy of a finite system is 
an analytic function and hence can be expanded about any point (e, m) in the parameter 
space. This is different in the Landau theory of the free energy of the infinite system. 
In this approximation the Helmholtz free energy is assumed to be analytic although a 
phase transition in the thermodynamic limit is related to a non-analytic point in the 
thermostatic potential. Consequently the critical behaviour obtained from the Landau 
approximation may differ from the true non-trivial behaviour of the physical systems. 
This is indeed the case below the upper critical dimension. 

The Taylor expansion of the microcanonical entropy of the finite system about the 
magnetisation m = (the equilibrium magnetisation of the high energy phase) yields 




(14) 
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the general series 

s{e,m) = s{e,m = 0) + ^ ap,,...,p„ie)Y[mf , (15) 

(pi,...,p„) i 

where the expansion coefficients are related to the derivatives of the entropy with respect 
to the magnetisations 

= ■ \ • • • d?n"Ae. m = 0) . (16) 

Pi - ■ ■ ■ Pn- 

The symmetry group Go of the entropy with respect to the magnetisation determines 
which coefficients ap^^,„^p^ appear in the expansion (fTH|l . Only the coefficients which 
correspond to monomials mFi . . . m^" that are invariant under the transformations 
g & Go are in general non-zero, all other coefficients necessarily vanish. In view of 
the symmetry Gq the expansion ()15|1 is in fact an expansion in terms of the symmetry- 
adapted harmonics of the group Gq. 

The physical behaviour for small magnetisations is already described by the low 
order terms in (fT3j) . Suppose the magnetisation is only one- dimensional and the 
microcanonical entropy is left invariant by the transformation m —m. Note that 
these are the symmetry properties of the entropy of the Ising model. Then only even 
powers contribute to the Taylor expansion and up to the fourth order term one gets 

s(e, m) = s(e,0) + -a(e)m^ + -b{e)m'^. (17) 

As the maximum of the entropy must not be at an infinite magnetisation the coefficient 
of the highest power of m has to be negative. This requirement is indeed fulffiled for 
an Ising model in the energy interval about Cc where a finite order parameter emerges. 
The Ising model is defined by the Hamiltonian 

n = -J2^^^^ (18) 

where the spins (Tj can be in the state +1 or —1 and the summation runs over neighbour 
pairs of spins. Figure ^ displays the derivatives dmis{e,0) of a two-dimensional Ising 
model with 200 spins near the energy Cc- The underlying two-dimensional lattice has 
a square-lattice topology. The density of states of the system is evaluated numerically 
exact by a microcanonical transfer matrix method (Binder, 1972; Creswick, 1995) so 
that it is possible to work out higher order derivatives reliably. The condition Q for 
thermostatic equilibrium leads to the equation 
d 

— — s(e, m) = m (a(e) + 6(e)m^) = . (19) 
am 

for msp(e). This relation has always the trivial solution msp(e) = and in addition the 
solutions 
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with a non-trivial, finite order parameter if tlie coefficient a(e) is positive. The stabihty 
condition 

^s(e,m,p(e))<0 (21) 

is always satisfied for the non-trivial solution (j^U)) provided it exists. The solution 
msp(e) = is only stable, i.e. corresponds to a maximum of the entropy, if the coefficient 
a(e) is negative. The coefficient a(e) changes its sign at the energy ec- Hence the 
equilibrium magnetisation of the entropy (fT7j) is given by 



msp(e) = ±e(a(e))^ 



a(e) 



b{e) 

where the Heaviside function is denoted by 0. As the coefficient 6(e) is negative in 
the vicinity of Cc and the two functions a and b are smooth the order parameter varies 
continuously as a function of the energy e. Below the energy Cc it attains a finite value. 
Note that at the energy Cc the solution m = becomes instable and the curvature of the 
entropy along the magnetisation changes its sign. The existence of two stable solutions 
below Cc signals the spontaneous breakdown of the global m — * ±m symmetry of the 
physical system. The depth 

A(e) = s(e,msp(e))-s(e,0) (23) 

of the convex dip along the magnetisation for fixed energy can also be expressed in terms 
of the expansion coefficients of the entropy surface. Near the critical point of the finite 
system one obtains the approximation 

A(e) ^ -0(«(e))^ • (24) 

The Landau expansion of the entropy surface around the equilibrium magnetisation 
of the high energy phase provides an explanation of the classical exponents describing 
the singular behaviour of the microcanonical physical quantities of the finite system. 
Introducing the reduced energy 

e := e — Cc (25) 

one can expand the functions a(e) and 6(e) about the critical energy Cc- 

a{e) = -Ae + age^ + . . . (26) 

and 

6(e) = -B + hie + 62^' + • • • (27) 

where the constants A and B are positive. Plugging these expansions into relations (j22|l 
and (j24|) one gets the leading behaviour 



m,^{e) = ±Q{-e)xl-^y\e\ (28) 
for the order parameter and 

A{e) = Q{-e)^e' (29) 
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for the depth of the convex dip in the hmit of small reduced energies e. Thus the variation 
of the spontaneous magnetisation in finite systems is described by the classical critical 
exponent (3 = 1/2. In addition one can deduce the leading behaviour of the curvature 
of the entropy along the magnetisation. For energies above the critical point, i.e. e > 
one has 



This leads to the critical exponent 7 = 1 characterising the divergent susceptibility 
in finite microcanonical systems. Note that the slope of the curvature dmmS at the 
spontaneous magnetisation above and below the critical point differ by a factor two. 
This is also observed in the curves obtained from the numerically calculated entropy 
surface (see discussion below and figure El). From the definition (jH)) of the magnetic 
field and the expansion ()17|) it is apparent that the magnetic field at the critical energy 
Cc vanishes with an exponent 5 = 3 in the limit of small magnetisations. Similarly 
the critical exponent a can be shown to be for both definitions of the specific heat 
(compare and (fT^ ). 

The expansion of the entropy surface hence describes the asymptotic behaviour of 
the physical quantities in the vicinity of the critical point of the finite system. It has to 
be stressed once again that the expansion of the entropy surface is not an approximation 
like the Landau expansion of the Helmholtz free energy. The Taylor expansion of the 
microcanonical entropy about an arbitrary point in the (e, m) space is always possible 
due to the analyticity of the entropy of finite systems. For the infinite system such an 
expansion about the critical point can not be performed as the thermostatic potential is 
singular precisely at this point. Therefore non-trivial critical exponents can only emerge 
in the thermodynamic limit. 

Although the physical quantities of finite microcanonical systems are always 
characterised by classical exponents precursors of the non-trivial exponent of the infinite 
system show up if the evolution of the quantities is considered for a series of different 
system sizes. This evolution of the physical quantities with increasing size L can be 
investigated with effective critical exponents (Hiiller/Pleimling, 2002) or microcanonical 
finite size scaling relations (Kastner et al., 2000). The microcanonical finite size scaling 
theory is conceptionally similar to the corresponding finite size scaling theory of the 
canonical ensemble. 

At this stage one delicate point has to be stressed again. Whereas the Taylor 
expansion of the entropy can be carried out without any ambiguity for systems with 
continuous energies and magnetisations (e.g. xy model) this seems to be doubtful for 
discrete models such as the Ising model. A suitable chosen continuous function may 
be fitted to the data. This leads, however, to a certain degree of arbitrariness which is 
still related to the choice of this function. Alternatively, differentials can be replaced by 
(centre) differences leading to a discrete set of data points. Nevertheless, any properly 




(30) 



and for negative reduced energies one has 




(31) 
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chosen function will reproduce the change of the curvature along the magnetisation 
at m = and Cc, leading to maxima at non-zero magnetisations. It is precisely this 
bifurcation at Cc that gives rise to classical exponents. Different fit-functions will only 
lead to slightly different coefficients, the overall dependence of the expansion an e and 
m will be the same. The considerations in this section strictly valid only for systems 
with continuous e and m give a heuristic understanding of the classical behaviour of 
physical quantities of discrete systems. 

4. Entropy surface with Cnv symmetry 

Landau expansion of entropies with Cnv invariance 

A general method that allows the determination of the invariant homogeneous 
polynomials which can appear in the Taylor expansion of the entropy is based on the 
concept of the complete rational basis of invariants of a group. The determination of 
the Landau expansion of the free energy of an infinite system with the help of the 
complete rational basis of invariants was first proposed by Gufan in (Gufan, 1971) (see 
also Toledano/Toledano, 1987). In this subsection this general method is applied to a 
microcanonical entropy surface with a Cnv symmetry. A prominent example of a physical 
system with this symmetry property of the entropy surface is the n-state vector Potts 
model which will be investigated in section 14.21 for a two-dimensional lattice. 

An invariant homogeneous monomial of degree p with respect to a finite group can 
be constructed as a linear combination of products of a limited number of polynomials 
forming the complete rational basis of invariants. In polar coordinates mi = r cos and 
m2 = r sin the complete rational basis of invariants of the finite group Cnv consists of 
the two polynomials and cos ncj). The polynomial = + m\ is always invariant 
under the modulus preserving transformations of the group Cnv 

The entropy of a finite microcanonical system with the symmetry 



can be expanded into a Taylor series. In polar coordinates the polynomials appearing in 
the expansion are sums of products of the monomials of the rational basis of invariant 
of Cnv The basic monomials contain the angular dependence of the magnetisation 
exclusively in the form cos ncf). Thus the entropy in polar coordinates is of the form 



This form of the entropy already determines properties that are independent of the 
details of the Taylor expansion, as e.g. the highest degree in the truncated Taylor 
series. These properties are solely a consequence of the symmetry of the entropy. With 
h = cosncf) the equilibrium equations are 



s(e,mi,m2) = s(e, C„^(mi, m2)) 



(32) 




(33) 



d_ 

dr 



s(e, r,b) = 



(34) 
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for the modulus of the spontaneous magnetisation and 

d ^ 

sin(n0)— s(e, r, 6) = 



(35) 



for its angular dependence. The equilibrium condition for the angular part can be 
satisfied in two ways. Either the derivative di,s vanishes or sinn0 is zero. The latter is 
the case for (p = kn/n with k = 0,...,2n — 1. The former way to satisfy equation ()35j] 
is discussed in section 14.31 The stability of the solution of (jHlj) and (jH^ requires that 
the corresponding Hessian is a negative definite matrix. Suppose that the stability is 
guaranteed for the radial order parameter component, i.e. drr-s < 0. At = kir/n the 
Hessian in polar coordinates reduces to 



and hence the negativity of H is satisfied if {—!)'' dbS > 0. The factor (—1)^ can either 
be +1 or —1 leading to maxima or saddle points of the microcanonical entropy surface 
at the solutions of the equilibrium conditions (jMj) and For a positive dbS the 

spontaneous magnetisation can have the orientations = kn/n with even k. In this 
situation the extrema on directions with odd k are saddle points of the entropy surface. 

These considerations show that the direction of the spontaneous magnetisation of 
an entropy with Cnv invariance is already determined from the symmetry properties with 
respect to the magnetisation. The extrema of the entropy appear along the symmetry 
lines of the associated regular n-polygon in the magnetisation plane. 

To get further information about the behaviour of the entropy surface for small 
magnetisations a Taylor expansion must actually be performed. The precise form of 
this Taylor expansion, i.e. the expansion coefficients and the degree of the truncated 
Taylor polynomial determines the stability of the extrema and the form of the energy 
dependence of the actual order parameter rsp(e) = |msp(e)|. To ensure the stability in 
angular direction one must have dbS ^ and hence the Taylor expansion has to contain 
at least the nth degree monomial r'^fc = r"^ cos n^. 

Entropy of the four-state vector Potts model 

The general results of the previous subsection can be illustrated for an entropy with 
symmetry in the magnetisation space. The rational basis of invariants is then given by 
{r^, cos 40} and the Taylor expansion of the entropy s{e, r, 0) up to the fourth degree 
term has the general form 



To ensure the appearance of a stable extremum for finite magnetisations the coefficients 
of the fourth degree term have to satisfy B{e) + C(e) cos 40 < for all angles 0. The 
extrema of the entropy are along the symmetry directions = kir/A with k = 0, . . . ,7 




(36) 




(37) 
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of the square in the two-dimensional plane. The stability in angular direction requires 
that the inequality 

C{e) cos k7r>0 (38) 

is satisfied. For a positive coefficient C{e) this has the consequence that the maxima of 
the microcanonical entropy surface lie on the directions G {0, 7r/2, tt, 37r/2} provided 
the stability of the radial order parameter is ensured. The extrema on the four remaining 
symmetry directions correspond to saddle points in the entropy function. 

An example of a classical spin model with C^^ symmetry is the vector Potts model 
(Wu, 1982) with four states whose Hamiltonian is given by 

n = -J2cos(^^ih-k,)'^ . (39) 

(id) 

The spins ki can take on the values 1, 2, 3 and 4 and be visualised by unit vectors 
with the orientations 0, it/2, it and 37r/2 in a two-dimensional plane. The intensive 
magnetisation of the system with spins is given by 

mi = cos (-ki) (40) 



N ^ \2 

i 
i 

The system has four equivalent ground states with the spontaneous magnetisations 
(1,0), (0,1), (—1,0) and (0,-1). These four ground state magnetisations define a 
square in the magnetisation plane. 

For a two-dimensional system with a finite number of spins the entropy has one 
single maximum at zero magnetisation above the critical energy Cc. Below Cc the 
entropy has four equivalent maxima along the directions defined by the coordinate 
axes. The extrema on the diagonals of the coordinate system correspond to saddle 
points of the entropy surface. In figure El the level curves of the entropy of a two- 
dimensional system with 64 spins are shown for two energies below the critical one. The 
data are calculated with the transition observable method that allows a highly efficient 
and accurate determination of microcanonical entropy surfaces (Hiiller/Pleimling, 2002). 
The findings show that the angular equilibrium equation ()35j) is satisfied by the condition 
sin 40 = for the four-state vector Potts model. The appearance of stable phases along 
the coordinate axes suggests that the coefficient C{e) is indeed positive. However it is 
hardly possible to obtain a reliable estimate of the fourth derivative from an entropy that 
has been calculated by means of Monte Carlo simulations. A direct determination of 
the expansion coefficients B{e) and C{e) is therefore not possible. To get an impression 
of the values of these coefficients one can use an indirect approach. The coefficient 
A[e) of the second degree term in the expansion (P7|l is determined by differentiating 
the entropy. The coefficients -B(e) ± C(e) corresponding to the direction and 7r/2 are 
then varied to achieve a good description of the data by the fourth degree expansion for 
small magnetisations. For the energy —0.75 near Cc = —0.7 this gives for the system 
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with L = 8 the value —1.14 for B and the positive value 0.04 for C in agreement with 
the general considerations above. For a system with L = 10 this procedure yields the 
coefficients B = —0.804 and the positive value 0.048 for C at the energy —0.86 below 
the critical energy —0.806 ± 0.005. Note however that this procedure may lead to large 
errors due to the subjective assessment of the agreement of the Taylor polynomial with 
the data for small magnetisations. The simulated entropy for the system with 100 spins 
together with its fourth degree Taylor approximation along the mi direction is displayed 
in figure El 

The behaviour of the spontaneous magnetisation of the two-dimensional system 
with 64 spins resulting from the variation of the position of the maximum is displayed 
in figure El Its energy dependence near the critical energy Cc = —0.7 is most suitably 
characterised by a square root behaviour, i.e. by the exponent f3 = 1/2. This can 
be understood from the general considerations in section El The non-zero spontaneous 
magnetisation msp{e) defines a vector in the order parameter space with a energy- 
dependent modulus |m5p(e)| and a fixed, energy-independent direction mo (compare 
relation Once the direction tuq has been chosen, e.g. along the mi axis, the 

entropy depends only on the modulus |m| leading to an entropy function with a one- 
dimensional order parameter. This function is just the cut of the entropy surface along 
the mo direction. The expansion of this function with respect to |m| can be performed 
as in section El and yields the characteristic asymptotic square root energy behaviour 
of the microcanonical equilibrium magnetisation of finite systems near the transition 
energy Cc. The curvature parallel to the spontaneous magnetisation at the equilibrium 
macrostate of the system with 100 spins near the critical point is shown in figure It 
results in a diverging susceptibility characterised by the critical exponent 7 = 1. The 
amplitudes from the left and from the right differ approximately by a factor 2. These 
observations are in agreement with the general findings obtained by Taylor expanding 
the microcanonical entropy surface. Note however that the susceptibility of the four- 
state vector Potts model is a tensor. Due to the symmetry ()32j) of the entropy surface 
the susceptibility tensor is diagonal with a non-vanishing parallel and perpendicular 
component with respect to the equilibrium magnetisation. The behaviour of the depth 
of the convex dip of the entropy surface below the critical energy Cc is shown in figure El 
for the system with 100 spins. Once again the parabolic energy dependence suggested 
by the above considerations is confirmed. 

4-3. Energy- dependent angular magnetisation 

To complete the discussion of the angular equilibrium equation this section focuses 
on the possible solution 9^5 = 0. The derivative d^s contains the radial magnetisation 
so that the two equilibrium equations for rgp and 0sp are coupled. This set of equations 
can be solved for h and gives h^p or equivalently 0sp as a function of r^p and hence it will 
depend on the energy. Consequently The direction of the spontaneous magnetisation 
exhibits an energy dependence. This dependence on the energy is infiuenced by the 
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precise form of the truncated Taylor expansion which determines the order parameter 
as a function of the energy. The stabihty of the extremum requires additionally that the 
expansion must contain monomials up to the degree 2n at least. To see this consider 
the component H^^ of the Hessian which is given by 

Hj^i = !^{smn<P,^fd,,s (42) 

at the spontaneous magnetisation. This requirements for the Taylor expansion has the 
further consequence that the phase diagram with respect to the expansion coefficients 
of the entropy allows only an isolated continuous transition (Gufan/Sakhnenko, 1973). 

The author does not know of any model system that exhibits an energy-dependent 
orientation of the microcanonically defined order parameter. Nevertheless, the general 
statements of this section can be illustrated with a model entropy which is invariant 
under the symmetry group 6*4^. With b = cos 40 and e = e — Cc the model entropy to 
investigate is given by 

sis, r, b) = soie) - ^er^ - ^r" + ^e^r^ - ^b\^ (43) 
2 4 4 8 

with positive coefficients A, B and G below Cc. The coefficient C has to be non-zero. 

For the following discussion it is assumed to be negative. The equilibrium equations for 

a non-zero spontaneous magnetisation (rgp, 0sp) are 

-Ae- Br^ + Cs^br'^ - Gb'^r^ = (44) 

and 

sin 40 (CeV^ - Gbr^) = . (45) 
The second equation can be solved for b which yields 

At this stage it is already obvious that the direction of the magnetisation vector in 
equilibrium will exhibit a dependence on the energy. Plugging this result into (jUj) gives 
rise to the expression 



|kl (47, 

for the radial order parameter. In view of equation ()46p one obtains 

GB^ 

fesp = cos40,p = -^^e . (48) 

Inverting this result for an angle 0sp G [0,7r/4] one ends up with the energy-dependent 
spontaneous direction 

1 

4™iGA 

Note that a different interval for the angle has to be chosen for a positive coefficient G. 
The resulting curve (rsp(e), 0sp(£:)) in the two-dimensional mim2-plane is schematically 



sp(£:) = - arccos ( -p^e ] . (49) 
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depicted in figure [7| Only one solution for the angle of the spontaneous magnetisation 
has been considered so far. As the entropy surface (pSj) is invariant under the group 
C^y the other solutions can be obtained by applying the transformations of C^v onto the 
solution (rsp, 0sp) from above. This results in a star of eight possible order parameters 
with an energy-dependent direction. The appearance of eight solutions is related to the 
eight different choices of the interval of length tt for the values of 40 in the interval 
[0, Svr] when inverting relation (PHj) . The Hessian of the model entropy evaluated at the 
solution for the order parameter is diagonal. The stability of the solution requires that 
the this Hessian has negative diagonal entries. This is indeed the 

= -2A\e\ - (50) 

Gr 

and 

Hf^ = -4{sm4<P^,yGr% (51) 
are both negative. 

5. Summary 

The appearance of classical critical exponents in physical observables of finite 
microcanonical systems such as the order parameter and the susceptibility is related 
to the assumed analyticity of the microcanonical entropy of finite systems. This can be 
demonstrated by Taylor expanding the entropy function about the well-defined critical 
point of the finite microcanonical system. The general findings are confirmed by concrete 
investigations of spin systems. This suggests in turn that the assumption of an analytic 
microcanonical entropy function is indeed justified although a proof is still lacking. An 
interesting example of a classical spin system is the four-state vector Potts model with 
a two-component order parameter. The form of the Taylor expansion of the entropy 
of the four-state vector Potts model is determined by the order parameter symmetry. 
The general consequences of the C^y symmetry on the form of the expansion of the 
microcanonical entropy and the resulting behaviour of the physical quantities for small 
magnetisations can be nicely verified by Monte Carlo calculations. The analyticity 
together with the knowledge of the order parameter symmetry allows already far- 
reaching statements about the physics of finite systems in the microcanonical ensemble. 
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Figure 1. The derivatives d^2s{e, 0) and d„^is{e, 0) of the entropy surface of a 
two-dimensional Ising model with 200 spins in a square lattice geometry. The 
fourth order derivative is multiplied by a factor / = 0.1. The second order 
derivative changes its sign at the transition energy Cc = —0.934 whereas the 
fourth order derivative stays negative in the vicinity of Cc- 




Figure 2. Level curves of the density of states of the four-state vector Potts 
model with 64 spins in two dimensions for the two energies —0.91 (a) and —1.06 
(b) below the critical point Cc = —0.7 it 0.005. Four equivalent maxima appear 
along the bisectors of the square defined by the ground state magnetisations. 
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Figure 3. The entropy of the two-dimensional four-state vector Potts model 
with 100 spins at the energy —0.86. The filled circles display the simulated data 
along the mi direction. The Taylor approximation 1)3 7|) up to the fourth degree 
term is shown as a dashed line. The second degree coefficient 0.0257 is obtained 
by differentiation, the fourth degree coefficient j{B+C) = —0.189 by variation. 
Additional higher order terms are necessary to describe the behaviour of the 
entropy for large magnetisations. 




Figure 4. The spontaneous magnetisation of the four-state vector Potts model 
with 64 spins in two dimensions as a function of the energy. The data suggest 
a critical energy Cc = —0.7 ± 0.005. 
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Figure 5. The (equilibrium) curvature of the entropy surface of the two- 
dimensional vector Potts model with 100 spins along the mi direction. Near 
the critical point Cc = —0.806 ± 0.005 the curvature vanishes leading to a 
divergent susceptibility. At the critical point the moduli of the gradients from 
the left and from the right differ by the factor 1.98 in good agreement with the 
prediction from the Taylor expansion. 



<^ 0.006 




Figure 6. The depth of the convex dip of entropy surface of the two- 
dimensional vector Potts model with 100 spins along the mi direction. Near 
the critical point Cc = —0.806 it 0.005 the energy variation is most suitably 
described by a parabola. 
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Figure 7. Schematic behaviour of the magnetisation curve (sohd hne) (rgp, (psp) 
in the two-dimensional plane that results from the model entropy H43|) . The 
angles —ir/S and vr/S are indicated by dashed lines. A second solution 
(dotted line) for the order parameter is obtained by reflecting the spontaneous 
magnetisation about the mi axis. Application of the rotations in 6*4^ to these 
two solutions gives the remaining solutions for the order parameter. 



